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Collagen is an important extracellular matrix component that 
directs many fundamental cellular processes including differen- 
tiation, proliferation and motility. The signalling networks driving 
these processes are propagated by collagen receptors such as 
the f31 integrins and the DDRs (discoidin domain receptors). 
To gain an insight into the molecular mechanisms of collagen 
receptor signalling, we have performed a quantitative analysis of 
the phosphorylation networks downstream of collagen activation 
of integrins and DDR2. Temporal analysis over seven time points 
identified 424 phosphorylated proteins. Distinct DDR2 tyrosine 
phosphorylation sites displayed unique temporal activation 
profiles in agreement with in vitro kinase data. Multiple clustering 
analysis of the phosphoproteomic data revealed several DDR2 
candidate downstream signalling nodes, including SHP-2 (Src 
homology 2 domain-containing protein tyrosine phosphatase 2), 
NCK1 (non-catalytic region of tyrosine kinase adaptor protein 1), 
LYN, SHIP-2 [SH2 (Src homology 2)-domain-containing inositol 



phosphatase 2], PIK3C2A (phosphatidylinositol-4-phosphate 3- 
kinase, catalytic subunit type 2a) and PLCL2 (phospholipase 
C-like 2). Biochemical validation showed that SHP-2 tyrosine 
phosphorylation is dependent on DDR2 kinase activity. Targeted 
proteomic profiling of a panel of lung SCC (squamous cell 
carcinoma) DDR2 mutants demonstrated that SHP-2 is tyrosine- 
phosphorylated by the L63V and G505S mutants. In contrast, 
the 163 8F kinase domain mutant exhibited diminished DDR2 and 
SHP-2 tyrosine phosphorylation levels which have an inverse 
relationship with clonogenic potential. Taken together, the results 
of the present study indicate that SHP-2 is a key signalling node 
downstream of the DDR2 receptor which may have therapeutic 
implications in a subset of DDR2 mutations recently uncovered 
in genome- wide lung SCC sequencing screens. 

Key words: cell signalling, collagen, discoidin domain receptor, 
lung cancer, mass spectrometry, phosphoproteomics. 
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INTRODUCTION 

Collagen is the most abundant protein in mammals and plays a 
critical role in maintaining tissue structural integrity [1]. There is 
an increasing appreciation that collagen has additional functions 
in propagating signal transduction networks which drive cellular 
behaviour in both physiological and pathological conditions. 
A number of different transmembrane collagen receptors have 
been identified, with the integrin family of receptors being 
the most well-studied [2]. Integrins are heterodimeric type 1 
transmembrane molecules composed of an a and subunit that 
engage a wide variety of extracellular matrix components and 
initiate intracellular signalling networks through the recruitment 
of downstream effector proteins [3,4]. The primary collagen- 
binding integrins are the a 1/3 1 and a 2/31 heterodimers, and 
these receptors have been shown to bind specific consensus 
sequence motifs in triple-helical regions on fibrillar collagen 
molecules [2,5]. In addition to the integrins, the DDRs (discoidin 
domain receptors) are another class of collagen-binding receptors. 
They differ from the integrins in that they belong to the RTK 



(receptor tyrosine kinase) superfamily and have intrinsic phospho- 
transferase activity through their cytoplasmic KD (kinase domain) 
[6-8]. Although classical growth factor-activated RTKs, such as 
the EGFR (epidermal growth factor receptor), display rapid and 
transient receptor activation dynamics [9], the DDRs are unique in 
that they exhibit delayed and sustained receptor phosphorylation 
upon binding to collagen [6,7]. The use of synthetic peptide 
toolkits to probe for receptor-binding specificity indicates that 
the DDRs bind to specific amino acid consensus sequences on 
fibrillar collagen that are distinct from integrin-binding motifs 
[10-12]. 

DDR2 is one of two members of the DDR family. It is highly 
expressed in mesenchymal cells and is important for a variety 
of developmental processes, in particular bone and cartilage 
formation [13]. DDR2 binds exclusively to fibrillar collagen and 
collagen X [14]. DDR2 also contributes to disease progression, 
including hepatic fibrosis, osteoarthritis and cancer [15-18]. 
Large-scale phosphoproteomic screens of human tumours have 
identified DDR2 to be highly phosphorylated in a subset of 
non- small-cell lung cancers, cholangiocarcinomas and sarcomas 
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Abbreviations used: CDK1, cyclin-dependent kinase 1; DDR, discoidin domain receptor; DMEM, Dulbecco's modified Eagle's medium; DYRK1A, 
dual-specificity tyrosine-phosphorylation-regulated kinase 1A; EGFR, epidermal growth factor receptor; ERK, extracellular-signal-regulated kinase; EV, 
empty vector; GO, Gene Ontology; HEK, human embryonic kidney; HRP, horseradish peroxidase; IL, interleukin; IMAC, immobilized metal-ion-affinity 
chromatography; KD, kinase domain; mAb, monoclonal antibody; MCAM, multiple clustering analysis methodology; NCK1 , non-catalytic region of tyrosine 
kinase adaptor protein 1 ; PIK3C2A, phosphatidylinositol-4-phosphate 3-kinase, catalytic subunit type 2a; PLCL2, phospholipase C-like 2; RFB, radiometric 
filter binding; RTK, receptor tyrosine kinase; SCC, squamous cell carcinoma; SHIP-2, SH2 (Src homology 2)-domain-containing inositol phosphatase 2; 
SHP-2, Src homology 2 domain-containing protein tyrosine phosphatase 2; SRM, selective reaction monitoring; TDA, template-directed assembly; TEAB, 
triethylammonium bicarbonate; TFA, trifluoroacetic acid. 
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[15,17,19]. Furthermore, cancer genome sequencing efforts have 
identified a series of oncogenic DDR2 point mutations that occur 
at low frequency in lung SCC (squamous cell carcinoma) [20,21]. 
A recent RNAi screen of protein tyrosine kinases has revealed that 
DDR2 is among a group of kinases that is required for extracellular 
matrix rigidity sensing and polarization in fibroblasts, implicating 
a role in mechanotransduction [22]. Although some of the 
functional roles of DDR2 in human health and disease have 
been identified, the signalling networks driving these biological 
processes remain largely unknown. 

In the present study, we have performed an unbiased 
quantitative MS-based phosphoproteomic analysis of collagen 
receptor signalling networks upon stimulation with collagen 
I. Temporal analysis over seven time points identified 424 
phosphorylated proteins. Using MCAM (multiple clustering 
analysis methodology), we find that a subset of phosphorylation 
sites on important signalling proteins such as SHP-2 (Src 
homology 2 domain-containing protein tyrosine phosphatase 2), 
NCK1 (non-catalytic region of tyrosine kinase adaptor protein 1), 
LYN, SHIP-2 [SH2 (Src homology 2)-domain-containing inositol 
phosphatase 2], PLCL2 (phospholipase C-like 2) and PIK3C2A 
(phosphatidylinositol-4-phosphate 3-kinase, catalytic subunit 
type 2a) strongly cluster with DDR2 phosphorylation dynamics, 
implicating these proteins as candidate downstream effectors of 
DDR2 signalling. A survey of lung SCC mutants of DDR2 using 
targeted proteomics show that SHP-2 is tyrosine-phosphorylated 
in a subset of these mutants. The present study demonstrates 
that SHP-2 is a key downstream component of DDR2 signalling 
and highlights the utility of phosphoproteomics in providing new 
insights into cell-matrix signalling. 



EXPERIMENTAL 

Cell culture, transfection and selection of cells 

HEK (human embryonic kidney)-293 cells were obtained from 
the A.T.C.C. and cultured in DMEM (Dulbecco's modified 
Eagle's medium) with 10 % FBS, 2 mM glutamine, 100 units/ml 
penicillin and lOOmg/ml streptomycin in an atmosphere of 
95% air/5% C0 2 at 37 °C. For expression of DDR2, HEK- 
293 cells were transfected with pcDNA3.1-DDR2 [11,23] or 
an empty vector control using the calcium phosphate method 
and selected in 400 /xg/ml zeocin (Invivogen) as described 
previously [24,25]. For DDR2 mutants, the QuikChange® XL 
site-directed mutagenesis kit (Stratagene) was used according to 
the manufacturer's protocol. Primers used for the mutagenesis 
are detailed in Supplementary Table SI (at http://www.biochemj. 
org/bj/454/bj4540501add.htm. Selected cells were pooled and 
analysed for DDR2 expression by immunoblotting analysis. 



Immunoblotting 

HEK-293-DDR2 or mutant cells were stimulated with 20 /xg/ml 
acetic-acid-solubilized rat tail collagen I (Sigma) at the indicated 
time points and lysed with RIPA lysis buffer [25 mM Tris/HCl 
(pH7.6), 150 mM NaCl, 1% Nonidet P40, 1% sodium 
deoxycholate and 0.1 % SDS] supplemented with protease and 
phosphatase inhibitors (Thermo Pierce) at 4°C. Equal amounts 
of protein, as determined by the bicinchroninic acid protein 
assay (Thermo Pierce), was applied in each lane of 4-12% 
or 10% Bis-Tris gels (Invitrogen). Following one-dimensional 
separation and transfer on to PVDF membrane, the membrane 
was incubated overnight at 4°C with 1:1000 dilutions of goat 
anti-DDR2 (R&D Systems), mouse anti-phosphotyrosine [4G10 



(Millipore) or PY100 (Cell Signaling Technology)] or rabbit anti- 
SHP2 (Santa Cruz Biotechnology) antibodies, a 1:500 dilution 
of a rabbit anti-(SHP2 p-Tyr 542 ) (Cell Signaling Technology) 
antibody or with a 1:5000 dilution of a mouse anti-a -tubulin 
(Sigma). After a 1 h incubation with secondary antibodies 
[HRP (horseradish peroxidase) -conjugated anti-goat IgG, anti- 
rabbit IgG or anti-mouse IgG at a 1:10000 dilution (Stratech 
Scientific)], immunoreactive bands were visualized by enhanced 
chemiluminescence (pico-L; Thermo Pierce) and the blots were 
exposed to X-ray XAR film (Kodak). 



Flow cytometry 

HEK-293-DDR2 cells were grown in six-well plates for 24 h. 
Cells were then dissociated with non-enzymatic cell dissociation 
solution (Sigma) and resuspended in PBS containing 1 % BSA. 
The cells were incubated for 30 min on ice with primary mAbs 
(monoclonal antibodies) at 10 /xg/ml in 100 /xl of PBS/BSA. Cells 
were then washed three times with PBS/BSA and incubated 
with FITC-conjugated goat anti-mouse IgG (F-9006, Sigma) 
for 30 min on ice. After three washes as above, the cells 
were resuspended in 2% formaldehyde in PBS. Data were 
subsequently collected on a BD LSRFortessa cell analyser using 
BD FACSDiva software 6.0 (BD Biosciences), and further 
analysed on FlowJo software 7.6.4 (Tree Star). Mouse anti- 
(integrin a I) mAb, clone FB12, was purchased from Millipore 
Chemicon and mouse anti-(integrin al) mAb, clone AK7, was 
from AbD Serotec. Mouse anti-(integrin /31), clone P5D2, was 
purified from hybridoma cells obtained from the Developmental 
Studies Hybridoma Bank, University of Iowa, Iowa City, I A, 
U.S.A. 



Sample preparation and phosphopeptide immunoprecipitation 

HEK-293-DDR2 cells were washed with PBS and incubated 
overnight in serum- free media prior to stimulation with 20 /xg/ml 
acetic-acid-solubilized rat tail collagen I (Sigma) at the time points 
indicated. Cells were then lysed in 8 M urea, and subjected to 
reduction, alkylation and trypsin digestion as described previously 
[26]. Briefly, samples were reduced with 10 mM DTT for 1 h 
at 56 °C and alkylated with 55 mM iodoacetamide for 1 h at 
room temperature (20 °C). Samples were digested with 40 /xg 
of trypsin (Promega) overnight at room temperature. Peptides 
were desalted on a C i8 Sep-Pak Plus cartridge (Waters), eluted 
with 25 % acetonitrile and freeze-dried to dryness. Freeze-dried 
peptides were subjected to labelling with the iTRAQ 8-plex 
reagent (ABSciex) according to the manufacturer's instructions 
[27]. Enrichment of phosphotyrosine peptides was achieved using 
peptide immunoprecipitation as described previously [26,28]. 
Briefly, 30 /xl of Protein G Plus-agarose beads (Sigma) were 
incubated with 15 /xg of each of the anti-phosphotyrosine 
antibodies, pYlOO (Cell Signaling Technology) and 4G10 
(Millipore), in 200 /xl of immunoprecipitation buffer [100 mM 
Tris, 100 mM NaCl and 1% Nonidet P40 (pH7.4)] for 8h 
at 4°C. Beads were washed with rinse buffer [100 mM Tris 
and 100 mM NaCl (pH7.4)] and retained peptides were eluted 
from the antibody with 70 /xl of elution buffer [100 mM glycine 
(pH2.5)] for 30 min at room temperature. The eluate from 
the immunoprecipitation was subjected to IMAC (immobilized 
metal-ion- affinity chromatography) enrichment before MS 
analysis. After the depletion of phosphotyrosine-containing 
peptides, 200 /xg of the supernatant from the immunoprecipitation 
was subjected to two rounds of IMAC for global phosphorylation 
analysis. 
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IMAC phosphopeptide purification and LC-MS/MS analysis 

IMAC enrichment of phosphorylated peptides was performed 
using Ni-NTA (Ni 2+ -nitrilotriacetate) agarose beads (Qiagen) 
using a protocol adapted from [29]. A volume of 100 /xl (for 
phosphotyrosine peptides) or 200 jA (for global phosphopeptide 
enrichment) of beads were washed three times with water, and 
incubated with 100 mM EDTA (pH8.0) for 30min. The beads 
were then washed three times with water, and incubated with 
100 mM FeCl 3 for 45 min. After removing excess metal ions, 
beads were washed three times with water. The iTRAQ-labelled 
eluent from the immunoprecipitation was acidified with 25 jA 
of 10 % TFA (trifluoroacetic acid) and added to the iron-chelated 
resin. Then 200 jA of 80 % acetonitrile and 0.1% TFA was added 
to the mixture before incubation for 1 h at room temperature. 
Beads were washed twice with 80 % acetonitrile and 0.1 % TFA, 
followed by two washes with 80 % acetonitrile and 0.1% acetic 
acid, and two washes of 0. 1 % acetic acid. For the global phos- 
phorylation analysis, the flowthrough from the first IMAC was 
subjected to an additional round of IMAC enrichment. The beads 
from the three IMAC enrichments (one for the phosphotyrosine 
immunoprecipitation and two for the global enrichment) were 
separately eluted with either 40 jA (phosphotyrosine peptides) 
or 80 jA (global phosphopeptide enrichment) of 250 mM sodium 
phosphate (pH 8.0) for 30 min and loaded directly on to a reverse- 
phase (Ci 8 ) pre-column (100 jim internal diameter, packed with 
10 cm of 10 jim C i8 beads). The pre-column was attached to an 
analytical column (50 jim internal diameter fused silica capillary 
packed with 10 cm of 5 jim Q 8 beads) with an integrated 
electro spray bottleneck tip with an approximate 1 jim orifice. 
For the phosphotyrosine analysis, peptides were eluted using 
a 125-min gradient with solvents A (1 % acetic acid) and B 
(water/acetonitrile/acetic acid at 10:89:1, by vol.): 10 min from 
0 to 13% B, 95 min from 13 to 42% B, 10 min from 42 
to 60% B, and 10 min from 60 to 100% B. For the global 
phosphopeptide analysis, peptides were eluted using a 220-min 
gradient with: 10 min from 0 to 13 % B, 190 min from 13 to 42 % 
B, 10 min from 42 to 60 % B, and 10 min from 60 to 100 % B. 
Eluted peptides were directly electrosprayed into a QqTOF mass 
spectrometer (QSTAR Elite; ABSciex) operated in information- 
dependent acquisition mode. MS/MS spectra of the five most 
intense peaks with two to five positive charges in the MS scan were 
automatically acquired with previously selected peaks excluded 
for 60 (phosphotyrosine) or 90 (global phosphorylation) s. 



Data analysis 

MS/MS spectra were searched against a Homo sapiens protein 
database (NCBI) by using MASCOT (version 2.2; Matrix 
Science) with trypsin as the enzyme and allowing up to three 
missed cleavages. Oxidation of methionine and phosphorylation 
of serine, threonine, tyrosine were included as variable 
modifications (0.15 Da MS/MS tolerance and 2.2 Da peptide 
tolerance), while carbamidomethylation of cysteine and iTRAQ 
modification of the -NH 2 lysine side chain and the N-terminus 
were included as fixed modifications. Peptide sequence validation 
was further confirmed manually for each of the peptides identified 
by checking the raw MS/MS data for possible mixed spectra, 
non- as signed abundant peaks and phosphorylation position. 
Phosphopeptide quantification was determined via Protein Pilot 
(ABSciex) by calculating the peak area for iTRAQ marker ions. 
The Protein Pilot software corrects for isotopic contamination 
associated with iTRAQ reagents as the signal for each isotopic 
tag contributes to the signal of the other tags. Quantification 
results were additionally manually validated. Each condition was 



normalized against the 12 1 . 1 channel to obtain fold changes across 
all seven conditions. To account for protein loading differences in 
the seven samples, a small fraction (~0.1 %) of the supernatant 
from the tyrosine phosphopeptide immunoprecipitation was 
analysed by LC-MS/MS, thereby providing quantification for 
the non-phosphorylated peptides in each sample. Protein loading 
quantification was then used to normalize the iTRAQ marker ion 
data for phosphorylated peptides. 



SRM (selective reaction monitoring) 

For SRM assays, cell lysates were prepared as detailed above 
for iTRAQ experiments. Analyses were performed using the 
equivalent of the same amount of cell lysate (1.5-2.4 mg 
depending on experiment) per condition. Following digestion 
and Sep-Pak desalting, phosphotyro sine-containing peptides 
were immunoprecipitated using 10 jig of the pYlOO antibody 
and 30 jA of Protein G Plus-agarose beads (Calbiochem). 
Immunoprecipitated peptides were eluted in 40 jA of elution 
buffer [100 mM glycine (pH 2.5)] and beads were removed 
by centrifugation at 5000 g for 3 min. Eluted peptides were 
then transferred to a fresh tube, and 2 jA of a heavy peptide 
standard mix was added per sample to allow for normalization 
of precipitated endogenous peptide levels between runs. Heavy 
peptides sequences are detailed in Supplementary Table S2 (at 
http://www.biochemj.org/bj/454/bj4540501add.htm). 

Samples were analysed using a Q-Trap 4000 instrument 
(ABSciex). Samples containing heavy peptide standards were 
loaded on to a reverse-phase (C i8 ) pre-column (100 jim internal 
diameter, packed with 5-10 cm of 10 jim C i8 beads). The pre- 
column was attached to an analytical column (50 jim internal 
diameter fused silica capillary packed with 10 cm of 5 jim 
Ci 8 beads) with an integrated electrospray bottleneck tip with 
an approximate 1 jim orifice. Peptides were eluted using a 
75-min gradient with solvent A (1 % acetic acid) and B 
(water/acetonitrile/acetic acid at 10:89:1, by vol.): 10 min from 
0 to 10% B, 45 min from 10 to 34% B, 10 min from 34 to 
47% B, and 10 min from 47 to 100% B. Transitions were 
monitored for endogenous and heavy phosphopeptides as detailed 
in Supplementary Table S3 (at http://www.biochemj.org/bj/454/ 
bj4540501add.htm). Peak areas were calculated using Analyst 
vl.5 software, and peak areas for endogenous peptide transitions 
were normalized to the corresponding transitions for the spiked 
heavy peptide standards. 



Dataset preparation and MCAM implementation 

The PTMScout [30] interface was used to analyse the phospho- 
proteomic dataset. The parameters of clustering used in the initial 
MCAM round and those pruned from the final round are provided 
in Supplementary Table S4 (at http://www.biochemj.org/bj/454/ 
bj4540501add.htm). Set pruning was determined by removing 
those sets whose removal improved the overall biological 
enrichment by 10% or greater, while not decreasing the impact 
of any one biological category by more than 5 % . Statistical 
significance of biological enrichment was calculated using the 
PTMScout MCAM interface for GO (Gene Ontology) terms, 
Pfam domains, Scansite binding and kinase predictions, linear 
amino acid sequence motifs, and enrichment for the domain 
a phosphorylation site falls in. This analysis used the default 
PTMScout parameters of 5x10" 2 for an a-corrected value, 
corrected by the Hochberg and Benjamini FDR procedure [31], 
1 x 10 " 5 for significant Pfam domain consideration, 1 x 10 ~ 2 cut- 
off for motif search analysis, and a Scansite cut-off of 3. At the 
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time statistical significance of enrichment was calculated, GO 
terms on PTMScout were at version 1.2 and downloaded June 
26, 2012 from GO (http://www.geneontology.org). The Matlab 
code for enrichment analysis and pruning was downloaded from 
http://ptmscout.mit.edu/MCAM_July_2011.zip. Further details 
regarding the process of MCAM, including enrichment analysis 
and pruning, can be found in [32]. 



In vitro kinase assay 

Using a TDA (template- directed assembly) approach which 
organizes functional complexes of membrane-associated proteins, 
His-tagged DDR2 and/or Src (see Supplementary Figure 
SIC at http://www.biochemj.org/bj/454/bj4540501add.htm for 
sequences) were mixed at the indicated amounts with a 1:5000 
ratio of TD A/enzyme. Enzyme solutions (2x; 1.5 /xM) were 
prepared in assay buffer 1 [AB1; 20 mM Hepes (pH 7.2), 10 mM 
MgCl 2 and 0.01 % Tween 20] and mixed with equal volumes 
of either (i) 200 /xM ATP and 200 /xM Axltide peptide substrate 
(KKSRGDYMTMQIG) in AB 1 for loading into PAGE gels, or (ii) 
200 /xM ATP (with 1.5 /xCi of [y- 32 P]ATP) and 200 /xM peptide 
substrate in AB 1 for RFB (radiometric filter binding) assays. Each 
activity assay was incubated for 30min at room temperature. 
RFB assay samples were stopped by spotting 20 /xl on to P- 
81 phosphocellulose filters pre-wetted with 0.75% phosphoric 
acid. Filters were washed four times with 35 ml of 0.75 % 
phosphoric acid, then immersed in 3 ml of scintillant (Ultima 
gold; PerkinElmer) and counted on a Beckman LSC-1800 liquid 
scintillation counter. 

Gel assay samples were stopped by the addition of 
reducing/denaturing loading buffer, then boiled for lOmin, 
quickly cooled, and loaded into 30 /xl wells on a 10% precast 
Tris-glycine gel (Bio-Rad Laboratories). Gels were run, stained 
with Coomassie Brilliant Blue and destained. Excised gel bands 
were diced and destained in consecutive 50 % acetonitrile/50 mM 
TEAB (triethylammonium bicarbonate) and 50 mM TEAB 
washes until all of the stain was removed. The gel pieces were 
shrunk amd then reduced with lOmM DTT and alkylated with 
55 mM iodoacetamide. The proteins were subjected to digestion 
with modified porcine trypsin (Promega) and incubated at 37 °C 
overnight. After digestion, the proteolytic peptides were extracted 
with 50 % acetonitrile/formic acid and TEAB washes, the extract 
solution was reduced in volume to approximately 10 /xl in a speed- 
vac and subjected to LC-MS/MS sequencing analysis. 



ELISA 

Cell lines were serum-starved overnight before stimulation with 
collagen I (Sigma) or acid control at the time points indicated 
and lysed. For P5D2-blocking experiments, HEK-293-DDR2 
cells were pre-treated with P5D2 (10 /xg/ml, R&D Systems) 
for 30 min before stimulation with collagen for 24 h and lysis. 
Lysis buffer for extraction of DDR2, PLCL2, SHIP2 and SHP- 
2 was 1 % Nonidet P40, 20 mM Tris (pH 8.0), 137 mM NaCl, 
10% glycerol and 4.5 mM EDTA; for LYN lysis buffer was 
1 mM EDTA, 0.5 % Triton X-100 and 6 M urea in PBS (pH 7.4); 
and for NCK1 lysis buffer was 50 mM Hepes, 0.1 mM EDTA, 
0.1 mM EGTA, 120 mM NaCl and 0.5 % Nonidet P40 (pH 7.5). 
For determination of phospho-SHP2 Tyr 542 levels, cells were 
lysed in 50 mM Hepes, 5 mM EDTA, 0.1 mM EGTA, 0.5% 
Nonidet P40 and 120 mM NaCl (pH 7.5). All lysis buffers were 
supplemented with protease and phosphatase inhibitors (Thermo 
Pierce). An ELISA for phospho-DDR2 was performed using the 
Human Phospho-DDR2 DuoSet IC ELISA kit (R&D Systems). 



Phospho-SHP2 was determined using the Phospho-SHP2 (Y542) 
DuoSet IC ELISA kit (R&D Systems). All of the remaining 
phosphoproteins were measured by custom sandwich ELISAs 
using the following capture antibodies diluted in PBS: total SHP- 
2 capture antibody (2 /xg/ml) from the SHP-2 DuoSet IC kit (R&D 
Systems); anti-PLCL2 mAb clone 2D 10 (5 /xg/ml, Sigma); anti- 
LYN polyclonal antibody (5 /xg/ml, catalogue number ab77400, 
Abeam); anti-SHIP2 mAb, clone T. 194.8 (1 :500 dilution, Thermo 
Scientific); and anti-NCKl mAb, clone 714506 (3 /xg/ml, R&D 
Systems). ELISA plates (Costar) were coated with capture 
antibody overnight at room temperature and subsequently blocked 
with 1 % BSA in PBS for 2 h at room temperature. 

Each sample (100 /xl) was loaded in duplicate and incubated 
for 2 h at room temperature, followed by the addition of 100 /xl of 
HRP- conjugated phosphotyrosine mouse mAb (1:1000 dilution, 
pYlOO, Cell Signaling Technology) in diluent [20 mM Tris, 
137 mM NaCl, 0.05% Tween 20 and 0.1% BSA (pH7.4)] 
(2 h at room temperature). The reaction was visualized by the 
addition of 100 /xl of chromogenic substrate [TMB (3,3 / ,5,5 / - 
tetramethylbenzidine), R&D Systems] for 30 min and the reaction 
was stopped with 50 /xl of 1 M H 2 S0 4 . Absorbance at 450 nm was 
measured using an ELISA plate reader. Plates were washed five 
times with wash buffer [PBS (pH7.4) containing 0.1 % Tween 
20] after each step. 



Colony formation assay 

Cell lines were resuspended in a solution composed of 1 mg/ml 
rat tail collagen (BD Biosciences) in DMEM with the pH adjusted 
to 7.4 with 1 M NaOH. Cells (400) were seeded in each well of a 
96- well plate in 100 /xl of the collagen solution. After the collagen 
gel solidified, 100 /xl of DMEM (10 % FBS) was added to the top 
of the gel and media was changed every other day. Cells were 
grown for 15 days and colonies were counted. 



RESULTS 

Generation of cell line and MS experimental strategy 

To identify the phosphorylation-mediated signalling networks 
downstream of DDR2 activation, HEK-293 cells were transduced 
to stably express DDR2 (HEK-293-DDR2 cells) which were then 
pooled before immunoblotting to assay for DDR2 expression. The 
HEK-293 cell line has been extensively used as a model system to 
investigate DDR signalling [6,11,33]. Stimulation of HEK-293- 
DDR2 cells with acid-soluble collagen I resulted in the activation 
of DDR2 with delayed and sustained tyrosine phosphorylation 
over 24 h (Figure 1A) as reported previously [6]. These cells 
maintained high expression levels of the major collagen integrins, 
/3l, al and a 2 integrin subunits, as demonstrated by flow 
cytometry assessment of surface expression levels of these 
receptors (Figure IB). Consistent with a recent characterization 
of integrin levels in HEK-293 cells [34], HEK-293-DDR2 cells 
maintained high levels of ^1 integrin subunit. Additionally, the 
cells expressed significant levels of the a 2 subunit and lower levels 
of the a 1 subunit. 

In order to minimize any confounding signals that may be 
associated with growth factors commonly found in serum in cell 
culture media, the HEK-293 -DDR2 cells were serum- starved 
before stimulation with 20 /xg/ml collagen I over seven time 
points (0-24 h) (Figure IC). Cells were lysed and the samples 
were prepared in two biological replicates as described in the 
Experimental section. Peptides from the seven conditions were 
stable isotope-labelled with the iTRAQ isobaric reagent, mixed 
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Figure 1 Characterization of cell line and MS experimental strategy 

(A) Immunoblot of DDR2 activation over 24 h after stimulation with 20 /xg/ml collagen I. Phosphotyrosine (pY) is measured by the 4G10 antibody. (B) Integrin expression in HEK-293-DDR2 cells. 
The cells were stained on ice with 10 /xg/ml of the primary mAbs followed by FITC-conjugated goat anti-mouse IgG and analysis by flow cytometry. Dotted line, secondary antibody only; grey 
line, anti-(integrin a1); black line, anti-(integrin al)\ filled grey histogram, anti-(integrin £1). Shown are representative data for three experiments. (C) Experimental strategy for phosphoproteomic 
analysis. HEK-293-DDR2 cells were stimulated across seven time points after the addition of collagen I and lysed. Proteins were extracted, digested and isotopically labelled with the iTRAQ 8-plex 
reagent. Phosphotyrosine enrichment was performed by anti-phosphotyrosine immunoprecipitation followed by IMAC. Global phosphorylation was assessed by two rounds of IMAC enrichment of 
the supernatant (SN) of the immunoprecipitation. Eluted peptides were then analysed by LC-MS/MS and the resulting data were processed for phosphopeptide identification and quantification as 
described in the Experimental section. 



and tyrosine-phosphorylated peptides were immunoprecipitated 
with a mixture of pan-specific anti-phosphotyrosine antibodies. 
Following immunoprecipitation, the phosphotyrosine peptides 
were further enriched using IMAC and analysed by LC- 
MS/MS. After the sample was depleted of phosphotyrosine- 
containing peptides, a global enrichment of phosphopeptides 
was performed with the iTRAQ- labelled supernatant from the 
immunoprecipitation being subjected to two rounds of IMAC 
before LC-MS/MS. In total, quantitative phosphorylation profiles 
were identified for 646 phosphosites on 649 phosphopeptides 
on 424 proteins across the seven time points (see the list of 
protein phosphorylation sites at http://www.biochemj.org/bj/454/ 
bj4540501add.htm). Of the phosphopeptides that were identified, 
manual sequencing of mass spectra unambiguously localized the 
site of phosphorylation on 549 of the peptides. 



DDR2 receptor phosphorylation upon collagen I activation 

The phosphoproteomic analysis revealed six phosphorylated 
peptides containing five phosphosites on DDR2 that were 



dynamically regulated upon collagen activation (Figure 2A). Two 
of these sites (Tyr 736 and Tyr 740 ) are located on the activation 
loop of the receptor which contains three tyrosine residues 
(Tyr 736 /Tyr 740 /Tyr 741 ) that are highly conserved among a number 
of RTKs, including the insulin receptor and c-Met [26,35]. These 
two sites occurred either as a singly or a doubly phosphorylated 
peptide (Tyr 736 + Tyr 740 ) and displayed a sharp increase in 
phosphorylation levels from 8 to 24 h, which is in agreement 
with the delayed activation kinetics of this receptor [6]. Two 
additional novel phosphosites which have not been previously 
described and are found in the KD of the receptor, Tyr 684 and 
Tyr 813 , exhibited similar kinetics, with maximal phosphorylation 
achieved at 24 h. In contrast, the Tyr 481 site present on the 
juxtamembrane domain of the receptor displayed constitutive 
phosphorylation levels and had a distinct profile compared with 
the other receptor phosphorylation sites. 

To confirm that the conserved tyrosine sites on the activation 
loop are associated with DDR2 kinase activity, we performed 
an in vitro kinase assay using TDA of His-tagged DDR2 in 
the presence of Src [36]. It has previously been shown that 
Src is required for activation of DDR2 in an in vitro setting 
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Figure 2 DDR2 receptor phosphorylation and in vitro kinase assay 

(A) Plots of relative DDR2 phosphorylation levels as a function of time after collagen stimulation. Temporal profiles of DDR2 sites show distinct responses to collagen stimulation. Measurements 
are expressed relative to the 24 h time point. (B) In vitro kinase assay measuring the incorporation of 32 P into the Axltide substrate peptide. Equal amounts (750 nM) of Src and DDR2 were mixed 
together with 300 pM TDA in kinase assay buffer. Data represent the fold change relative to control (DDR2). (C) Immunoblot of total phosphotyrosine levels (using the 4G10 antibody) obtained from 
the in vitro kinase reaction showing increased phosphorylation levels that correlate with DDR2 kinase activity. 



Table 1 List of phosphorylation sites identified in the in vitro kinase assay 



Phosphosite 


Sequence 


Condition 










DDR2 


DDR2 + TDA 


DDR2 + Src 


DDR2 + Src + TDA 


DDR2 negative control 


Ser 446 


RMLDDEMTVSLpSLPSDSSMFNNNR 


+ 


+ 


+ 


+ 




Ser 461 


SSpSPSEQGSNSTYDR 


+ 


+ 


+ 


+ 


+ 


Tyr 471 


SSSPSEQGSNSTpYDR 






+ 


+ 


+ 


Tyr 481 


IFPLRPDpYQEPSR 


+ 


+ 


+ 


+ 




Ser 674 


HEPPNpSSSDVR 








+ 




Ser 696 


FMATQIApSGMK 


+ 


+ 






+ 


Tyr 740 


NLYSGDpYYR 








+ 




Tyr 741 


NLYSGDYpYR 








+ 





and when co-expressed in baculo virus [37]. The TDA approach 
has been used to restore biological relevant functionality to 
recombinant membrane-associated proteins, or fragments thereof 
[38-40]. TDA of DDR2 in combination with Src enhanced both 
the kinase activity as well as the total tyrosine phosphorylation 
levels of the receptor compared with Src or TDA treatment alone 
(Figures 2B and 2C). Src itself does not phosphory late the Axltide 
substrate used in these experiments (Supplementary Figure SI A). 
The resulting kinase assay reaction was run on a gel and the 
band corresponding to DDR2 (Supplementary Figure SIB) was 
excised and subjected to LC-MS/MS analysis to identify the 
phosphorylation sites on DDR2 that are present upon kinase 
activation (Table 1). The sequencing analysis confirmed that two 
sites (Tyr 740 and Tyr 741 ) on the activation loop were phosphory lated 
only upon the addition of Src and TDA at which maximal 
kinase activity was observed. In agreement with the cellular 
phosphoproteomic data, Tyr 481 was seen to be phosphory lated 
under all experimental conditions, even in its basal inactive state, 
indicating constitutive phosphorylation of this site. Interestingly, 
a number of previously undescribed serine phosphorylation sites 
(Ser 446 , Ser 674 and Ser 696 ) were also identified in this in vitro 
analysis, and additional future work will be required to determine 
the functional relevance of these phosphosites in cells. 



Multiple clustering analysis identifies components of the DDR2 
signalling network 

Since our phosphoproteomic analysis represents the activation 
of both DDR2 and the integrins, we subjected the data to 



MCAM to isolate the DDR2-specific signalling components. 
MCAM is an unsupervised learning approach that requires no 
previous assumptions of the network properties under study 
and applies a variety of clustering algorithms, number of 
target clusters (K), distance metrics and data transformations, 
in a combinatorial fashion. Following application of clustering, 
feature selection is employed in the form of statistically enriched 
biological annotations within clusters [32]. The unambiguous 
phosphorylation dataset was chosen for this analysis since the 
exact sites of phosphorylation were defined. After clustering 
and parameter refinement, 216 distinct clustering sets were 
obtained to create an ensemble of clustering sets. To identify 
candidate downstream components specific to DDR2 signalling, 
co-occurrence analysis was performed to establish the frequency 
at which pairs of phosphosites co-cluster robustly together in 
the ensemble of clustering sets. These data are represented by 
the co-occurrence matrix (Figure 3A). This matrix was then 
clustered to find groups of robustly co-clustered phosphoproteins 
to produce a final hard partition based on the ensemble results. 
This analysis was done using Hierarchical clustering and Ward 
linkage to produce 40 clusters. We found that the first seven 
clusters have a higher than average co-occurrence frequency and 
that the majority of DDR2 phosphosites are found in clusters 
1-3 in the co-occurrence matrix (Figure 3B, and Supplementary 
Table S5 at http://www.biochemj.org/bj/454/bj4540501add.htm). 
These phosphorylation sites included DDR2 Tyr 740 in cluster 
1, Tyr 736 , Tyr 736 + Tyr 740 and Tyr 684 in cluster 2, and Tyr 813 in 
cluster 3. Phosphory lated proteins in these clusters exhibit low 
phosphorylation levels at the early time points that reached a 
maximum at 24 h after collagen stimulation (Figure 3C). MCAM 
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Figure 3 MCAM analysis of the phosphoproteomic dataset 

(A) Co-occurrence matrix for the phosphoproteomic dataset. Co-occurrences are calculated by counting the number of times any two phosphopeptides cluster together across all 216 clustering 
sets. The matrix is then normalized as a percentage of the number of times they cluster out of 216 times and subsequently clustered using hierarchical clustering with Ward linkage. The first seven 
clusters that display higher-than-average co-occurrence are highlighted by the red box. (B) The co-occurrence map of the first seven clusters expanded from (A). Cluster 1 contains DDR2 Tyr 740 
as well as SHP2 Tyr 62 . Cluster 2 contains three DDR2 phosphopeptides containing Tyr 736 , Tyr 736 /Tyr 740 and Tyr 618 . Cluster 3 contains DDR2 Tyr 813 . (C) The dynamics of the members of each of 
the three DDR2 phosphorylation site-containing robust clusters. (D) Distribution of co-occurrences in the ensemble clustering result. A single hard clustering from the ensemble is obtained by 
cutting the Ward-linked hierarchically clustered co-occurrence matrix in (A) into 40 clusters. The co-occurrences within every cluster is then plotted to give an idea of the distribution. Comparison of 
phosphotyrosine/phosphotyrosine (pY/pY) sites are in red, phosphoserine/phosphothreonine (pS/pT):phosphoserine/phosphothreonine (pS/pT) are in blue, and cyan is of the co-occurrence events 
between phosphotyrosine (pY) and phosphoserine/phosphothreonine (pS/pT) sites. 



analysis highlights that despite these DDR2 phosphorylation sites 
having relatively similar dynamics, particular features of the 
individual sites can be distinguished. For example, the Tyr 740 
cluster is probably distinguished, by clustering, from the Tyr 813 
site by lower phosphorylation at early time points. Importantly, 
although there is a large degree of overlap between the DDR2- 
containing clusters (clusters 1-3), these clusters are distinct from 



the remainder of the dataset. Consistent with the role of DDR2 as 
an RTK, a plot of the individual pairwise co-occurrence values for 
the members in each cluster reveals that DDR2 phosphorylation 
sites robustly co-cluster with other tyrosine phosphosites in the 
dataset (Figure 3D). 

Robust co-clustering of DDR2 receptor phosphorylation with 
a large number of tyrosine-phosphorylated proteins in the 
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Figure 4 Visualization of the fold change in phosphorylation levels of important signal transduction adaptors and effector proteins upon temporal activation 
of DDR2 as determined by MCAM analysis 

Representative proteins were chosen from the first three clusters identified to co-cluster with high frequency with DDR2 receptor phosphorylation in the co-occurrence matrix analysis. 



dataset suggests that these proteins are likely to be candidate 
downstream components specific to the DDR2 signalling network. 
For instance, NCK1 is an adaptor protein that was found in 
cluster 2 and has previously been associated with DDR signalling 
[41]. Other members in these three clusters include important 
integrators of signal transduction, such as LYN, PIK3C2A, 
SHP-2, SHIP-2, ERK1 (extracellular- signal-regulated kinase 1) 
and PLCL2 (summarized in Figure 4). Acid control treatment 
of HEK-293-DDR2 cells did not lead to phosphorylation of 
DDR2 and its downstream components SHP-2, NCK1, LYN, 
PLCL2 and SHIP-2 as measured by ELIS A (Supplementary Fig- 
ure S2 at http://www.biochemj.org/bj/454/bj4540501add.htm), 
indicating that the observed increases in downstream protein 
phosphorylation are the result of collagen addition. 

Collagen stimulation results in the activation of both the 
integrins and DDR2. The activation of these two distinct classes 
of receptors may occur on different timescales. To establish that 
these candidate proteins are specific to DDR2 and not due to 
integrin activation, HEK-293-EV (empty vector) control cells 
which express endogenous levels of collagen-binding integrins, 
but not DDR2 [25], were subjected to collagen treatment at both 
early (30 min and 2 h) and late (24 h) time points. As shown in 
Supplementary Figure S2, these cells fail to display increases in 
the phosphorylation of SHP-2, NCK1, LYN and PLCL2 at all 
time points, providing support that integrin activation is not the 
driver of the observed increases in tyrosine phosphorylation. The 
exception is SHIP-2, which has a small phosphorylation increase 
at 24 h in the HEK-293-EV cells, but a greater enhancement 
in the DDR2-expressing cells. Finally, to ascertain if activation 
of both DDR2 and integrins is required for signalling network 
propagation, HEK-293-DDR2 cells were pre-treated with a 



/31-integrin-blocking antibody (P5D2) before stimulation with 
collagen for 24 h (Supplementary Figure S3 at http://www. 
biochemj.org/bj/454/bj4540501add.htm). Blocking integrin ac- 
tivation does not alter the tyrosine phosphorylation status of 
downstream signalling compared with control. These experiments 
demonstrate that the observed up-regulation of protein tyrosine 
phosphorylation identified by the MCAM analysis occurs 
independently of integrin activation and is the result of DDR2 
activation by collagen. 



SHP-2 tyrosine phosphorylation is dependent on DDR2 kinase 
activity 

Phosphorylation of SHP-2 at Tyr 62 , a protein tyrosine phosphatase 
that has previously been shown to be an important regulator of 
growth factor RTK signalling, was found to co-cluster robustly 
with DDR2 Tyr 740 (90%), and with Tyr 684 and Tyr 736 (70%) 
(see Supplementary Table S6 at http://www.biochemj.org/bj/454/ 
bj4540501add.htm for co-clustering frequency). We independ- 
ently confirmed this finding by performing SRM analysis of this 
phosphorylation site on SHP-2 upon stimulation of HEK-293- 
DDR2 cells with collagen. By spiking in a synthetic heavy phos- 
phopeptide of SHP-2 Tyr 62 in a targeted proteomics experiment, 
we show that SHP-2 is robustly tyrosine phosphorylated 24 h after 
collagen stimulation (Figure 5A). This observation was further 
validated using a phospho-specific antibody for SHP-2 Tyr 542 , a 
phosphorylation site that is important for normal ERK activation 
in response to growth factor signalling [42]. In agreement with 
the phosphoproteomic data, temporal stimulation of HEK-293- 
DDR2 cells with collagen I showed an increase in receptor 
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Figure 5 DDR2 phosphorylates SHP-2 in a temporal and kinase-dependent manner 

(A) SRM analysis of DDR2 activation loop (Tyr 736 and Tyr 740 /Tyr 741 ) and SHP-2 (Tyr 62 ) phosphorylation at 0 and 24 h post-stimulation with collagen I (n = 3). Values are means + S.E.M., **P < 0.01 
using paired Student's t test. Representative transitions are shown as indicated in Supplementary Table S3 (at http://www.biochemj.org/bj/454/bj4540501add.htm). (B) Immunoblot of DDR2 and 
SHP-2 (Tyr 542 ) phosphorylation shows a temporal up-regulation of SHP-2 phosphorylation upon simulation with collagen. Phosphotyrosine (pY) is measured by the 4G10 antibody. (C) Immunoblot 
of kinase-dead DDR2 mutants (K608M and K608E) indicate that tyrosine phosphorylation of SHP-2 (Tyr 542 ) is dependent on DDR2 kinase activity. EV, empty vector control. 



phosphorylation, accompanied by a concomitant increase in SHP- 
2 Tyr 542 phosphorylation levels (Figure 5B). Correlation analysis 
of the phosphorylation levels of Tyr 62 and Tyr 542 on SHP-2 across 
15 different measurements over a range of collagen stimulation 
time points by SRM and ELISA respectively (Supplementary 
Table S7 at http://www.biochemj.org/bj/454/bj4540501add.htm) 
demonstrates that there is a strong correlation between these 
two phosphorylation sites (Spearman correlation coefficient 
r = 0.9321, P < 0.0001). Importantly, the total abundance of the 
SHP-2 protein does not vary after 24 h exposure to collagen 
and thus protein abundance changes are not a contributing 
factor to the phosphorylation changes observed. SHP-2 tyrosine 
phosphorylation was dependent on DDR2 kinase activity, as 
HEK-293 cells expressing either control vector or a kinase- 
dead version of DDR2 (either K608M or K608E mutations) 
were unable to phosphorylate SHP-2 when cells were exposed 
to collagen I (Figure 5C). Taken together, these data indicate that 
SHP-2 is a downstream effector of DDR2 activation. 



SHP-2 is differentially tyrosine phosphorylated downstream of a 
subset of lung SCC-associated DDR2 mutants 

A number of cancer-associated DDR2 point mutations were 
recently identified at low frequency through cancer genome 
sequencing studies in lung SCC [21]. These mutations are 
found throughout the DDR2 gene and range from extracellular 



mutations to mutations in the KD [43]. The signalling pathways 
activated by these mutant receptors have not previously been 
characterized and we sought to determine whether SHP-2 was 
similarly phosphorylated downstream of these DDR2 mutants. 
HEK-293 cells stably expressing three mutants that are found 
in the extracellular collagen-binding discoidin domain (L63V), 
the juxtamembrane region (G505S) or the KD (163 8F) were 
engineered (Figure 6A). In the original study by Hammerman et 
al. [21], the authors demonstrated in NIH 3T3 cells and BaF3 cells 
that a subset of DDR2 mutants, including L63V and I638F, are 
oncogenic. However, these assays were performed in the absence 
of collagen and the authors did not provide any experimental 
evidence that DDR2 was activated and phosphorylated in the 
cell lines tested. Activation of DDR2 by fibrillar collagen 
has previously been shown to inhibit proliferation of human 
melanoma and fibrosarcoma cells in 2D and 3D assays [44]. To 
determine the phenotypic effect of DDR2 and its lung cancer 
mutants in HEK-293 cells in the context of its natural ligand, 
collagen, a colony formation assay measuring the clonogenicity 
of these cells in 3D collagen I gels was performed as described 
previously [45-47]. Consistent with the tumour- suppressive role 
of DDR2 in cancer cells induced by fibrillar collagen [44], 
we found that wild- type DDR2 expression greatly reduced 
colony formation compared with the control cells (Figure 6B). 
Interestingly, DDR2 mutants showed variable responses to 
collagen I exposure, with the L63V and G505S showing 
intermediate suppression, and the 163 8F mutant displaying no 
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Figure 6 Targeted proteomic profiling of lung SCC DDR2 mutants 

(A) Domain organization and location of DDR2 point mutations used in the present study. JM, juxtamembrane; TM, transmembrane domain. (B) Colony formation assay of DDR2 mutants grown in 
3D collagen I gels (n = 3). Values are means + S.E.M., statistical significance of wild-type and mutant DDR2 data compared with the empty vector control was performed by ANOVA with Dunnett's 
post-test where **P < 0.01 and *P < 0.05. (C) Immunoblot of mutant DDR2 cells after stimulation with collagen I for 24 h. Phosphotyrosine (pY) is measured by the pY100 antibody. (D) Normalized 
densitometry measurements of phosphotyrosine and phosphorylated SHP-2 blots (n = 3). ***P < 0.001 and **P < 0.01 . Normalization is performed relative to the loading control, tubulin. (E) SRM 
analysis of DDR2 receptor and SHP-2 phosphorylation in wild-type and lung SCC point mutations post activation with collagen I at 24 h [n = 3). Values are means + S.E.M., statistical significance 
of mutant DDR2 phosphorylation data compared with wild-type DDR2 was performed by paired Student's t test where ***P < 0.001 and *P < 0.05. Representative transitions are shown as indicated 
in Supplementary Table S3 at http://www.biochemj.org/bj/454/bj4540501add.htm. 
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significant effect on colony numbers compared with the control 
(Figure 6B). 

Western blot analysis of these cells confirmed that the L63V 
and G505S mutant receptors were tyrosine-phosphorylated in a 
similar fashion to wild-type DDR2, whereas the 163 8F mutant 
showed a significant reduction in receptor phosphorylation 
(Figures 6C and 6D). We further demonstrate that L63V and 
G505S were able to phosphorylate SHP-2 at Tyr 542 , but I638F was 
incapable of doing so. To characterize the differences in DDR2 
and SHP-2 tyrosine phosphorylation with site-specific resolution, 
wild-type and mutant cells were stimulated with collagen I and 
subjected to targeted SRM analysis. MS analysis revealed that 
distinct sites on DDR2 (Tyr 481 , Tyr 684 , Tyr 736 and Tyr 740 /Tyr 741 ) 
and SHP-2 (Tyr 62 ) displayed varied phosphorylation responses 
among the different mutants (Figure 6E). A common observation 
across all of the phosphorylation sites is that the 163 8F mutant 
significantly diminished tyrosine phosphorylation on DDR2 and 
SHP-2 in agreement with the Western blot analysis. 



DISCUSSION 

In the present study, we have performed an unbiased phospho- 
proteomic analysis of collagen receptor activation. Much of the 
published work on DDR2 has thus far focused on its extracellular 
domain and the nature of its binding to collagen [2,43]. The 
present study is, to our knowledge, the first comprehensive 
temporal study of collagen receptor signalling networks in cells. 
We have identified 424 phosphorylated proteins, including five 
tyrosine sites on the DDR2 receptor. Quantitative analysis of 
DDR2 receptor phosphorylation dynamics reveals that these sites 
are differentially regulated. The sites on the activation loop (Tyr 736 
and Tyr 740 ) were not phosphorylated in the absence of collagen 
and showed delayed phosphorylation dynamics when stimulated, 
culminating in a sharp increase in phosphorylation levels from 8 to 
24 h after collagen engagement. Also phosphorylated in a similar 
fashion are two novel receptor sites (Tyr 684 and Tyr 813 ) in the KD. 
Additionally, these sites were not identified in the in vitro kinase 
assay (Table 1), which suggests that they may be phosphorylated 
by kinases other than Src or DDR2 autophosphorylation. In 
contrast, Tyr 481 on the juxtamembrane domain of the receptor 
displays constitutive phosphorylation which was confirmed 
in vitro with recombinant DDR2 which was phosphorylated 
even in its basal state. Owing to the lack of specific reagents 
(e.g. antibodies) targeted against these phosphorylation sites on 
DDR2, the present study represents the first demonstration that 
individual sites on DDR2 are regulated in a distinct temporal 
fashion, which may have consequences on the recruitment of 
downstream effectors and ultimately cellular behaviour. 

MCAM analysis of the phosphoproteomic data identified 
a number of tyrosine-phosphorylated proteins that clustered 
robustly with DDR2 receptor phosphorylation (Supplementary 
Tables S5 and S6). Since collagen is known to stimulate an array 
of different receptors [2], this clustering approach was particularly 
useful for identifying proteins that are candidate downstream 
effectors of DDR2 signalling. Several of these effector proteins 
have been associated with DDR2 function. For instance, PI3K 
(phosphoinositide 3-kinase) and ERK activation is required for 
DDR2-dependent IL (interleukin)-6 secretion in primary human 
chondrocytes [48]. In addition, DDR2 increases the expression of 
matrix metalloproteinase-13 in synovial fibroblasts via the ERK 
pathway [49]. We focused our validation efforts on SHP-2, a 
protein tyrosine phosphatase that has previously been implicated 
in DDR1 signalling in MDCK (Madin-Darby canine kidney) 
cells [50]. This DDR1 study did not determine the tyrosine 



phosphorylation status of SHP-2 upon collagen stimulation. Using 
SRM and Western blot analysis, we demonstrate that SHP-2 
is tyrosine-phosphorylated at two sites (Tyr 62 and Tyr 542 ) upon 
DDR2 activation. We further show that DDR2 kinase activity is 
required for SHP-2 phosphorylation, providing additional support 
that SHP-2 is a specific downstream effector of DDR2. It should 
be noted that the clustering analysis that we have performed 
is correlative in nature and serves to highlight the candidate 
protein phosphorylation sites that are co-regulated with DDR2 
receptor phosphorylation. Although such clustering approaches 
are useful tools for hypothesis generation, whether these co- 
regulated events represent direct interactions between specific 
DDR2 phosphotyrosine sites and downstream candidate proteins 
requires additional experimental validation, such as the generation 
of tyrosine-to-phenylalanine mutants of DDR2 to probe for site- 
specific effects on SHP-2 phosphorylation. 

We have recently shown that the collagen-binding integrins are 
active in the HEK-293-DDR2 cells used in the present study [25]. 
However, as we did not observe any phosphorylation sites on 
the collagen-binding integrin subunits in the phosphoproteomic 
dataset, we are unable to definitively identify the signalling 
events attributed to integrin activation using MCAM. A previous 
global phosphoproteomic analysis of integrin-mediated adhesion 
of HeLa cells identified 517 phosphorylation sites that were 
regulated within 2h of adhesion on collagen-I-coated surfaces 
[51]. Two of these tyrosine phosphorylation sites [CDK1 
(cyclin-dependent kinase 1) Tyr 15 and DYRK1 A (dual-specificity 
tyrosine-phosphorylation-regulated kinase 1A) Tyr 321 ] were also 
found in our dataset, but displayed very different profiles from 
the integrin study. In the integrin study, CDK1 Tyr 15 showed 
increasing phosphorylation levels reaching a maximum at 2 h after 
collagen stimulation, whereas DYRK1A Tyr 321 displayed high 
basal phosphorylation levels that steadily decreased over time. 
Both of these sites showed constitutive phosphorylation levels 
that did not vary with time in our dataset. These differences could 
be the result of the distinct cell lines used in the two studies, the 
use of soluble collagen compared with collagen-coated surfaces, 
or the influence of DDR2 expression on integrin activation. 

Unlike lung adenocarcinomas, SCCs of the lung are poorly 
characterized and there are currently no effective targeted 
therapies to treat this disease. The identification of point mutations 
in the DDR2 gene in lung SCC by Hammerman et al. [21] 
suggests that this RTK may be a novel oncogenic target for 
therapy. However, their study was performed using classical 
transformation assays including anchorage dependence of NIH 
3T3 cells grown in soft agar and IL- 3 -independent growth of 
BaF3 cells. In contrast, results from the present study demonstrate 
that HEK-293 cells expressing DDR2 have reduced clonogenic 
potential compared with control cells when grown in 3D collagen 
gels. Inhibition of cancer cell growth by fibrillar collagen is 
well documented and our results are consistent with published 
data indicating that DDR2 inhibits cancer cell proliferation when 
challenged with fibrillar collagen [44,52,53]. Loss of DDR2 has 
also been shown to predispose hepatic tissues to colon carcinoma 
cell growth and metastasis [54]. Interestingly, the present study 
shows that the 163 8F KD mutant restores the growth potential 
of these cells to control levels. The tyrosine phosphorylation 
status of the DDR2 mutants has not previously been established. 
We employ targeted proteomics to determine, with site-specific 
resolution, the phosphorylation profiles of DDR2 and SHP-2 in 
mutant-expressing cells. We find that the L63V and G505S 
mutants are signalling- competent and are capable of driving both 
receptor and SHP-2 tyrosine phosphorylation. The 163 8F KD 
mutant conversely displays a reduced ability to mediate DDR2 
and SHP-2 phosphorylation. Taken together, our data suggests 
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that in the presence of its natural ligand, DDR2 acts as a growth 
suppressor, and one potential mechanism of action of the 163 8F 
KD mutant is to alleviate this suppression by eliminating receptor 
and SHP-2 phosphorylation (Figures 6B and 6D). However, the 
relationship between growth suppression and signal transduction 
is likely to be complex and multivariate in nature since the G505S 
mutant displays high receptor and SHP-2 phosphorylation, but 
only has a small effect on growth suppression. The contrasting 
oncogenic and tumour-suppressive properties of DDR2 are likely 
to be context-dependent. The lung is an important site of collagen 
deposition and modification (such as by cross-linking) during 
tumour progression and the tumour-suppressive properties of 
DDR2 may be more pertinent in this context [55,56]. 

In conclusion, the results of the present study highlight the use 
of phosphoproteomics in combination with clustering techniques 
to establish a high-resolution map of DDR2 signalling effectors. 
We further demonstrate that SHP-2 is a major downstream 
component of DDR2 signalling and is phosphorylated in a subset 
of DDR2 lung SCC mutants. These findings provide an insight 
into the nature of DDR2 networks and signalling modulations 
that occur upon acquisition of cancer-associated mutations. We 
have previously shown that distinct oncogenic mutations in 
EGFR deploy differential pathway utilization in cancer cells, 
which has therapeutic implications in the clinic [26,27,57]. The 
results of the present study provide further support for this 
concept where the 163 8F KD mutant down-regulates receptor and 
SHP-2 phosphorylation, and alleviates DDR2-mediated growth 
suppression. Six kinase doman mutations have been identified 
in lung SCC, and three additional mutations in endometrial 
carcinoma and colorectal cancer [21]. Our findings suggest that 
restoring DDR2 receptor activity and its downstream signalling 
effectors, such as SHP-2, may be important for preserving tumour 
suppression in a subset of cancer-associated KD mutants. 
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Figure S1 Details of the DDR2 in vitro kinase assay 

(A) In vitro kinase assay measuring the incorporation of 32 P into the Axltide substrate peptide. Src and DDR2 were mixed at a ratio of 1:20 together with TDA at a 1:10000 (TDA/enzyme) ratio in 
kinase assay buffer. Src does not phosphorylate the Axltide substrate peptide. (B) Coomassie Brilliant Blue-stained gel of in vitro kinase reaction samples that were subjected to LC-MS/MS analysis. 
(C) Sequence of recombinant DDR2 and Src that were used in the in vitro kinase assays. 
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Figure S2 ELISA measurements of tyrosine phosphorylation levels of DDR2 and its downstream effectors (SHP2, SHIP-2, PLCL2, LYN and NCK1) in 
HEK-293-DDR2 cells at 0, 0.5, 2 and 24 h post collagen I stimulation (n = 4) 

Values are means + S.E.M. with ***P < 0.001, **P < 0.01 and *P < 0.05, indicating a significant difference between t = 0 and t = 24 h as determined by paired Student's t test. As negative 
controls HEK-293-EV cells treated with collagen I and HEK-293-DDR2 cells treated with acetic acid were used. 
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Figure S3 ELISA measurements of tyrosine phosphorylation levels of DDR2 and its downstream effectors (SHP2, SHIP-2, PLCL2, LYN and NCK1) in 
HEK-293-DDR2 at 24 h post collagen I stimulation (n = 4) in the presence or absence (control) of the P5D2 integrin ySI -blocking antibody 

Values are means + S.E.M. There are no statistical differences between the P5D2 and control samples. 
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Table S1 Primer sequences for generation of DDR2 mutants 

The bases shown in bold indicate the nucleotide that was mutated from the corresponding 
wild-type to generate mutants using site-directed mutagenesis methodology. 



Table S2 Heavy phosphopeptide sequences used in the SRM assay 



Mutant 


Direction 


Primer sequence (5'— >3') 


L63V 


Forward 


CTGCCAAATATGGAAGGGTGGACTCAGAAGAAGGG 




Reverse 


CCCTTCTTCTGAGTCCACCCTTCCATATTTGGCAG 


G505S 


Forward 


GAGTCAGGCTGCAGCAGTGTTGTGAAGCCAG 




Reverse 


CTGGCTTCACAACACTGCTGCAGCCTGACTC 


K608M 


Forward 


GTCCTGGTGGCTGTGATGATGCTCCGAGCAGATG 




Reverse 


CATCTGCTCGGAGCATCATCACAGCCACCAGGAC 


K608E 


Forward 


TGTCCTGGTGGCTGTGGAGATGCTCCGAGCAGATG 




Reverse 


CATCTGCTCGGAGCATCTCCACAGCCACCAGGACA 


I638F 


Forward 


CTCAAGGACCCAAACATCTTCCATCTATTAGCTGTGT 




Reverse 


AC AC AG CTAATAG ATG G AAG ATGTTTG G GTC CTTG AG 



Phosphosite Peptide 



Amount spiked per sample 



DDR2 Tyr 481 IFPLRPD[Tyr(P0 3 H 2 )]QEPS[Arg( 13 C6; 15 N4)] 500 fmol 

DDR2 Tyr 684 TVS[Tyr(P0 3 H 2 )]TNL[Lys( 13 C6; 15 N 2 )] 100 fmol 

DDR2 Tyr 736 NL[Tyr(P0 3 H 2 )]SGDYY[Arg( 13 C 6 ; 15 N 4 )] 100 fmol 

DDR2 Tyr 740 NLYSGD[Tyr(P0 3 H 2 )]Y[Arg( 13 C 6 ; 15 N4)] 200 fmol 

SHP2 Tyr 62 IQNTGD[Tyr(P0 3 H 2 )]YDLYGGE[Lys( 13 C 6 ; 15 N 2 )] 2.5 pmol 



Table S3 Transitions used for the SRM assay 

Values in bold indicate representative transitions used in the Figures. 



Phosphosite 


Peptide 


Iransition (U1/U3) 


Fragment ion ID 


Collision i 


DDR2 Tyr 481 (heavy) 


IFPLRPD[Tyr(P0 3 H 2 )]QEPS[Arg( 13 C6; 15 N4)] 


569.9/369.2 


Y3 


36.6 






569.9/869.3 


Y6 


36.6 






569.9/498.3 


Y4 


36.6 


DDR2 Tyr 684 (heavy) 


TVS[Tyr(P0 3 H 2 )]TNL[Lys( 13 C 6 ; 15 N 2 )] 


507.2/813.3 


Y6 


28.7 






507.2/201.1 


B2 


28.7 






507.2/726.3 


Y5 


28.7 






507.2/483.3 


Y4 


28.7 


DDR2 Tyr 736 (heavy) 


NL[Tyr(P0 3 H 2 )]SGDYY[Arg( 13 C 6 ; 15 N 4 )] 


620.7/1013.3 


Y7 


35.0 






620.7/770.3 


Y6 


35.0 






620.7/511.3 


Y3 


35.0 






620.7/683.3 


Y5 


35.0 


DDR2 Tyr 740 (heavy) 


NLYSGD[Tyr(P0 3 H 2 )]Y[Arg( 13 C 6 ; 15 N 4 )] 


620.7/1013.3 


Y7 


35.0 






620.7/391.2 


B3 


35.0 






620.7/591.1 


Y3 


35.0 






620.7/706.2 


Y4 


35.0 


SHP2Tyr 62 (heavy) 


IQNTGD[Tyr(P0 3 H 2 )]YDLYGGE[Lys( 13 C 6 ; 15 N 2 )] 


912.4/789.4 


Y7 


53.7 






912.4/674.4 


Y6 


53.7 






912.4/561.3 


Y5 


53.7 






912.4/398.2 


Y4 


53.7 


DDR2 Tyr 481 (endogenous) 


IFPLRPD[Tyr(P0 3 H 2 )]QEPSR 


566.6/359.2 


Y3 


36.6 






566.6/859.3 


Y6 


36.6 






566.6/488.3 


Y4 


36.6 


DDR2 Tyr 684 (endogenous) 


TVS[Tyr(P0 3 H 2 )]TNLK 


503.2/805.3 


Y6 


28.7 






503.2/201.1 


B2 


28.7 






503.2/718.3 


Y5 


28.7 






503.2/475.3 


Y4 


28.7 


DDR2 Tyr 736 (endogenous) 


NL[Tyr(P0 3 H 2 )]SGDYYR 


615.7/1003.3 


Y7 


35.0 






615.7/760.3 


Y6 


35.0 






615.7/501.2 


Y3 


35.0 






615.7/673.3 


Y5 


35.0 


DDR2 Tyr 740 (endogenous) 


NLYSGD[Tyr(P0 3 H 2 )]YR 


615.7/1003.3 


Y7 


35.0 






615.7/391.2 


B3 


35.0 






615.7/581.1 


Y3 


35.0 






615.7/696.2 


Y4 


35.0 


SHP2 Tyr 62 (endogenous) 


IQNTGD[Tyr(P0 3 H 2 )]YDLYGGEK 


908.4/781.37 


Y7 


53.7 






908.4/666.4 


Y6 


53.7 






908.4/553.3 


Y5 


53.7 






908.4/390.2 


Y4 


53.7 
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Table S4 Dataset preparation and MCAM implementation 

Parameters used for the initial MCAM analysis and pruned from the final round of feature selection. FFT, fast Fourier transform. 



Parameter 



Original parameters of clustering 



Removed during feature selection 



K 5,8,11,14,17,20,23,26,29 5 

Transform Centre, FFT, differential, no transform, z-score, normMax, rangeScale, square root, pareto Centre, FFT, differential 

Distance Correlation, Euclidean, cityblock, cosine, Chebychev Correlation 

Algorithm Hierarchical, k-means, affinity propagation, self-organizing maps Hierarchical 

Number of cluster sets 950 21 6 (remaining) 



Table S5 Phosphorylation sites in the top seven clusters in MCAM analysis 

The sites in bold highlight the SHIP-2 and DDR2 phosphorylation sites in each cluster. 



Cluster 1 


Cluster 2 


Cluster 3 


Cluster 4 


Cluster 5 


Cluster 6 


Cluster 7 


PFDN6 Tyr 82 


SF3B1 4 Tyr 86 


BAIAP2 Tyr 491 


EDC4 Ser 723 


MYL9 Thr 19 


AAK1 Ser 637 


PRKAR1A Ser 83 


BANF1 Tyr 43 


DDR2 Tyr 684 


SCRIB Tyr 1360 


SFRS15 Ser 154 


PYGLSer 15 


ZC3H13 Ser 64 


TPR Ser 2155 


SHIP-2 Tyr 986 


SDCBP Tyr 50 


SNX9 Tyr 269 


N0C2LSer 672 ,Ser 673 


SHIP-2 Tyr 1135 


SRRM2 Ser 1987 


SRRM2 Ser 2100 , Thr 2104 


UTRNTyr 3111 


EPS15L1 Tyr 74 


TUBBTyr 340 


CCDC88ATyr 1799 


HSPA1 A Tyr 611 


SRRM2 Thr 1208 


TP53BP1 Ser 380 


MAGOHTyr 123 


DDR2 Tyr 736 


RNPS1 Tyr 205 


MAPK14 Tyr 182 


AKT1S1 Ser 183 


SRRM2 Ser 1694 


MAPI S Ser 657 


SHP2 Tyr 62 


ITSN1 Tyr 1132 


IGF2BP2 Tyr 40 


TBC1D15Ser 227 


TRIM28 Tyr 517 


RPLP1 Ser 101 , Ser 104 


SRRM2 Ser 2272 


PSAT1 Tyr 346 


NCK1 Tyr 112 


DDR1 Tyr 796 


CLPB Ser 23 


PTPRA Tyr 798 


KLC2 Ser 581 


TNKS1BP1 Ser 691 


LCP1 Tyr 28 


IGF2BP3 Tyr 39 


ERK1 Tyr 204 


RPLP1 Ser 101 


SYKTyr 323 


SPTAN1 Ser 1217 


SRRM2 Thr 1492 


PIN4Tyr 122 


PLCL2 Tyr 784 


VIM Tyr 117 


CRKRS Ser 685 , Ser 681 






STMN1 Ser 16 , Ser 25 


EDC4 Ser 729 


DDR2 Tyr 736 , Tyr 740 


LYN Tyr 306 


JUN Ser 73 






SRRM2 Ser 1179 


DDR2 Tyr 740 


AGFG1 Tyr 327 


HSPA1 A Tyr 41 


SPAG9 Thr 217 








CLTCTyr 1487 


ELM 02 Tyr 48 


LPPTyr 296 


SFRS1 Tyr 189 








SF3B14 Tyr 61 


ANKRD39 Tyr 65 


DDR2 Tyr 813 










FASN Tyr 130 


ACBD3 Ser 43 


PABPC1 Tyr 54 










ITSN1 Tyr 1054 


ADD1 Tyr 407 













PIK3C2ATyr 1595 
TRAP1 Tyr 498 
DDR1 Tyr 792 
GTF2E1 Tyr 91 
AK2 Tyr 190 
TUBGCP3 Tyr 114 
PSMC3 Tyr 132 
ACTBTyr 218 
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Table S6 Co-occurrence frequency listed by the DDR2 phosphorylation sites 

The sites in bold highlight SHP-2 Tyr 62 phosphorylation and their respective co-occurrence frequency with specific DDR2 phosphorylation sites. 



DDR2 Tyr 481 


DDR2 Tyr 684 




DDR2 Tyr 736 




DDR2 Tyr 740 




DDR2 Tyr 736 , Tyr 740 




DDR2 Tyr 813 




Phosphorylation Co-occurrence 


Phosphorylation 


Co-occurrence 


Phosphorylation 


Co-occurrence 


Phosphorylation 


Co-occurrence 


Phosphorylation 


Co-occurrence 


Phosphorylation 


Co-occurrence 


site frequency 


site 


frequency 


site 


frequency 


site 


frequency 


site 


frequency 


site 


frequency 


DDR2 Tyr 481 1.00 


DDR2 Tyr 684 


1.00 


DDR2 Tyr 736 


1.00 


PFDN6 Tyr 82 


1.00 


DDR2 Tyr 736 Tyr 740 


1.00 


DDR2 Tyr 813 


1.00 


CRKLTyr 207 0.67 


AGFG1 Tyr 327 


0.91 


SF3B1 4 Tyr 86 


0.94 


DDR2 Tyr 740 


1.00 


ITSN1 Tyr 1132 


0.88 


ERK1 Tyr 204 


0.84 


TPI1 Ser 21 0.62 


NCK1 Tyr 112 


0.90 


NCK1 Tyr 112 


0.91 


PIK3C2ATyr 1595 


0.98 


ELM02 Tyr 48 


0.88 


LYN Tyr 306 


0.82 


RANBP1 Ser 60 0.55 


DDR2 Tyr 736 


0.86 


IGF2BP3 Tyr 39 


0.91 


DDR1 Tyr 792 


0.97 


EPS15L1 Tyr 74 


0.83 


LPPTyr 296 


0.82 


COIL Ser 301 0.51 


IGF2BP3 Tyr 39 


0.86 


SDCBP Tyr 50 


0.87 


LCP1 Tyr 28 


0.97 


ANKRD39 Tyr 65 


0.83 


SCRIB Tyr 1360 


0.81 


WDR75 Ser 796 0.50 


SDCBP Tyr 50 


0.83 


DDR2 Tyr 684 


0.86 


GTF2E1 Tyr 91 


0.97 


RNPS1 Tyr 205 


0.80 


RNPS1 Tyr 205 


0.80 




ANKRD39 Tyr 65 


0.82 


AGFG1 Tyr 327 


0.86 


ACTBTyr 218 


0.97 


ADD1 Tyr 407 


0.79 


DDR1 Tyr 796 


0.80 




EPS15L1 Tyr 74 


0.82 


ADD1 Tyr 407 


0.85 


BANF1 Tyr 43 


0.95 


NCK1 Tyr 112 


0.78 


SNX9 Tyr 269 


0.78 




SF3B1 4 Tyr 86 


0.81 


ITSN1 Tyr 1054 


0.84 


EDC4 Ser 729 


0.93 


SCRIB Tyr 1360 


0.78 


TUBBTyr 340 


0.77 




ADD1 Tyr 407 


0.81 


ANKRD39 Tyr 65 


0.84 


CLTC Tyr 1487 


0.93 


IGF2BP3 Tyr 39 


0.77 


BAIAP2 Tyr 491 


0.73 




ELM 02 Tyr 48 


0.77 


EPS15L1 Tyr 74 


0.83 


SF3B14 Tyr 61 


0.92 


DDR1 Tyr 796 


0.77 


IGF2BP2 Tyr 40 


0.72 




ITSN1 Tyr 1132 


0.76 


ACBD3 Ser 43 


0.81 


MAGOH Tyr 123 


0.92 


ERK1 Tyr 204 


0.77 


VIM Tyr 117 


0.69 




ACBD3 Ser 43 


0.74 


ITSN1 Tyr 1132 


0.79 


UTRNTyr 3111 


0.90 


LYN Tyr 306 


0.77 


PABPC1 Tyr 54 


0.68 




PSAT1 Tyr 346 


0.74 


ELM 02 Tyr 48 


0.79 


SHP2 Tyr 62 


0.90 


TUBBTyr 340 


0.74 


DDR2 Tyr 736 , Tyr 740 


0.67 




ITSN1 Tyr 1054 


0.73 


PIN4Tyr 122 


0.78 


TRAP1 Tyr 498 


0.89 


SDCBP Tyr 50 


0.74 


ELM 02 Tyr 48 


0.64 




CLTCTyr 1487 


0.71 


PLCL2 Tyr 784 


0.76 


PSMC3Tyr 132 


0.88 


SNX9 Tyr 269 


0.73 


HSPA1 A Tyr 41 


0.63 




SHIP-2Tyr 986 


0.70 


PSAT1 Tyr 346 


0.73 


TUBGCP3 Tyr 114 


0.87 


LPPTyr 296 


0.72 


SDCBP Tyr 50 


0.61 




TUBGCP3 Tyr 114 


0.70 


AK2 Tyr 190 


0.73 


PSAT1 Tyr 346 


0.85 


DDR2 Tyr 736 


0.71 


ITSN1 Tyr 1132 


0.60 




SHP2 Tyr 62 


0.69 


TUBGCP3 Tyr 114 


0.72 


SHIP-2 Tyr 986 


0.85 


SF3B1 4 Tyr 86 


0.70 


ADD1 Tyr 407 


0.60 




DDR2 Tyr 736 , Tyr 740 


0.69 


DDR2 Tyr 736 , Tyr 740 


0.71 


AK2 Tyr 190 


0.77 


ACBD3 Ser 43 


0.70 


SF3B1 4 Tyr 86 


0.59 




TRAP1 Tyr 498 


0.69 


FASN Tyr 130 


0.70 


FASN Tyr 130 


0.75 


DDR2 Tyr 684 


0.69 


IGF2BP3 Tyr 39 


0.59 




DDR1 Tyr 792 


0.69 


SHIP-2Tyr 986 


0.70 


ITSN1 Tyr 1054 


0.75 


VIM Tyr 117 


0.67 


PLCL2 Tyr 784 


0.59 




AK2 Tyr 190 


0.69 


SHP2 Tyr 62 


0.69 


PIN4Tyr 122 


0.71 


DDR2 Tyr 813 


0.67 


EPS15L1 Tyr 74 


0.58 




PIN4 Tyr 122 


0.69 


TRAP1 Tyr 498 


0.69 


AGFG1 Tyr 327 


0.70 


IGF2BP2 Tyr 40 


0.65 


ANKRD39 Tyr 65 


0.58 




SF3B1 4 Tyr 61 


0.69 


CLTCTyr 1487 


0.68 


DDR2 Tyr 684 


0.68 


AGFG1 Tyr 327 


0.65 


DDR2 Tyr 736 


0.57 




FASN Tyr 130 


0.69 


SF3B1 4 Tyr 61 


0.65 


PLCL2 Tyr 784 


0.65 


PABPC1 Tyr 54 


0.63 


NCK1 Tyr 112 


0.57 




PIK3C2ATyr 1595 


0.68 


DDR1 Tyr 792 


0.65 


SF3B14 Tyr 86 


0.64 


PLCL2 Tyr 784 


0.63 


JUN Ser 73 


0.57 




GTF2E1 Tyr 91 


0.68 


GTF2E1 Tyr 91 


0.64 


DDR2 Tyr 736 


0.64 


HSPA1 A Tyr 41 


0.62 


ACBD3 Ser 43 


0.57 




PFDN6 Tyr 82 


0.68 


PFDN6 Tyr 82 


0.64 


ACBD3 Ser 43 


0.63 


ITSN1 Tyr 1054 


0.59 


PIN4Tyr 122 


0.55 




DDR2 Tyr 740 


0.68 


UTRNTyr 3111 


0.64 


IGF2BP3 Tyr 39 


0.58 


PIN4Tyr 122 


0.59 


ITSN1 Tyr 1054 


0.53 




ACTBTyr 218 


0.67 


DDR2 Tyr 740 


0.64 


SDCBP Tyr 50 


0.58 


BAIAP2 Tyr 491 


0.56 


SFRS1 Tyr 189 


0.53 




UTRNTyr 3111 


0.67 


PIK3C2ATyr 1595 


0.64 


NCK1 Tyr 112 


0.58 


AK2 Tyr 190 


0.50 


SFRS15 Ser 154 


0.52 




LCP1 Tyr 28 


0.66 


PSMC3Tyr 132 


0.63 


ANKRD39 Tyr 65 


0.55 






CCDC88ATyr 1799 


0.52 




BANF1 Tyr 43 


0.65 


ACTBTyr 218 


0.63 


EPS15L1 Tyr 74 


0.55 






SPAG9 Thr 217 


0.52 




PLCL2 Tyr 784 


0.65 


LYN Tyr 306 


0.62 


ADD1 Tyr 407 


0.55 






DDR2 Tyr 684 


0.50 




MAGOH Tyr 123 


0.65 


BANF1 Tyr 43 


0.62 


ITSN1 Tyr 1132 


0.51 












EDC4 Ser 729 


0.65 


MAGOH Tyr 123 


0.62 


ELM02 Tyr 48 


0.51 












PSMC3Tyr 132 


0.64 


LCP1 Tyr 28 


0.62 


JUN Ser 73 


0.51 












LYN Tyr 306 


0.58 


EDC4 Ser 729 


0.61 


SFRS1 Tyr 189 


0.50 












SNX9 Tyr 269 


0.57 


DDR2 Tyr 813 


0.57 
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I UDD lyl 
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RNPS1 Tyr 205 


0.55 


ERK1 Tyr 204 


0.56 
















DDR1 Tyr 796 


0.55 


DDR1 Tyr 796 


0.56 
















SCRIB Tyr 1360 


0.53 


RNPS1 Tyr 205 


0.55 
















ERK1 Tyr 204 


0.53 


SNX9 Tyr 269 


0.54 
















VIM Tyr 117 


0.50 


TUBBTyr 340 


0.52 
















DDR2 Tyr 813 


0.50 


LPPTyr 296 


0.50 
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Table S7 Correlation analysis for phosphorylation of SHP-2 Tyr 62 and Tyr 542 

HEK-293-DDR2 cells were treated with 20/xg/ml collagen I and harvested across a range of 
different time points. For each time point, equivalent lysates were harvested for both SRM and 
ELISA experiments. Each sample represents a time point after collagen stimulation of cells. Data 
have been normalized to sample 15. Spearman correlation coefficient r = 0.9321 , P < 0.0001 . 



Sample ELISA (Tyr 542 ) SRM (Tyr 62 ) 



1 


0.137 


0.063 


2 


0.362 


0.283 


3 


0.521 


0.403 


4 


0.541 


0.336 


5 


0.620 


0.360 


6 


0.655 


0.525 


7 


0.699 


0.391 


8 


0.809 


0.457 


9 


0.848 


0.618 


10 


0.969 


1.870 


11 


1.099 


0.664 


12 


1.107 


0.719 


13 


1.177 


2.891 


14 


1.868 


3.083 


15 


1.000 


1.000 
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